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Abstract. The linear stability of a shocked isothermal accretion flow onto a black hole is investigated in the 
inviscid limit. The outer shock solution, which was previously found to be stable with respect to axisymmetric 
perturbations, is, however, generally unstable to non-axisymmetric ones. Eigenmodes and growth rates are ob- 
tained by numerical integration of the linearized equations. The mechanism of this instability is based on the cycle 
of acoustic waves between their corotation radius and the shock. It is a form of the Papaloizou-Pringle instability, 
modified by advection and the presence of the shock. As such it can be generalized to non isothermal shocked 
accretion flows. Blobs and vortices are generated by the shock as a by-product of the instability. 
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1. Introduction 

Hydrodynamic instabilities of shocked accretion flows may 
explain some of the properties of X-ray binaries, such 
as their time variability. The structure of stationary ac- 
cretion flows involving shocks was described by Fukue 
(1987) and Chakarabarti (1989a, b). Even with the sim- 
ple inviscid hypothesis, the structure of shocked accre- 
tion flows is complex, and their stability is not yet fully 
understood. As noted by Nakayama (1992), the calcula- 
tions of Chakrabarti ( 1989a, b) is a study of the forced 
oscillations of the flow rather than an analysis of its in- 
trinsic stability. Nakayama (1992, 1993) introduced a new 
type of global instability between the inner sonic point 
and the shock. He found that, of the two possible shock 
positions, the inner one is unstable due to post-shock ac- 
celeration, while the outer one is stable due to post-shock 
deceleration. His conclusion was confirmed by Nobuta and 
Hanawa (1994), whose numerical simulations showed that 
the inner shock is completely destroyed by perturbations, 
while the outer one is stable. All the above works, how- 
ever, were based on axisymmetric calculations. Molteni, 
Toth & Kuznetsov (1999, hereafter MTK) performed 2- 
D simulations of an adiabatic flow with an outer shock 
and found a non-axisymmetric instability. They showed 
that the instability saturates at a low level, and a new 
asymmetric configuration develops, with a deformed shock 
rotating steadily. MTK pointed out that this effect may 
have relevant observational consequences, such as quasi- 
periodic oscillations (QPO). The mechanism of the insta- 
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bility was not explained by MTK, who briefly mentioned 
a possible link with the numerical simulations of Blaes & 
Hawley (1988). The Papaloizou-Pringle instability (1984, 
hereafter PPI) simulated by Blaes & Hawley is known to 
take place in discs or tori, in which the radial velocity is 
initially zero, whereas the flow simulated by MTK involves 
radial advection, an inner sonic point and an outer shock. 
The effect of advection on the PPI was investigated by 
Blaes (1987), who found that the PPI is strongly stabilized 
by advection at the inner boundary. The interpretation of 
the results of MTK in terms of the PPI is thus not obvious 
a priori. What is more, the flow studied by MTK is also 
potentially unstable by the advective-acoustic mechanism 
(Fogfizzo & Tagger 2000, Foglizzo 2001, 2002), based on 
the cycle of entropy /vorticity perturbations and acoustic 
waves in the subsonic region between a stationary shock 
and a sonic surface. 

The aim of this study is thus to understand the insta- 
bility mechanism at work in shocked accretion flows. For 
the sake of simplicity, the present linear stability anal- 
ysis is focused on isothermal accretion flows with con- 
stant angular momentum. The paper is organised as fol- 
lows. Linearized equations and boundary conditions are 
described in Section 2 and solved numerically in Section 
3. The instability mechanism is analysed in Section 4, and 
compared to the simulations by MTK in Section 5. 



2. Equations 

An inviscid, isothermal flow around a black hole is 
considered, in the pseudo-Newtonian potential intro- 
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Fig. 1. Radial profile of the Mach number in the unperturbed 
flow, for which I — 1.87 and Cs = 0.1. The two circles denote 
the inner and the outer sonic points, respectively. The dashed 
arrow shows the shock at rsh = 6.8. 

duccd by Paczynski and Wiita (1980), $ = -GM/{r - 
Tg). Equations are made dimensionless by using the 
Schwarzschild radius and the speed of Ught as reference 
units, i.e., = 1 and c = 1. In these units, the Keplerian 
frequency is noted Qk = ^/{r — l)(2r)5. In order to 
overcome the technical difficulty of treating a realistic 
vertical structure, three types of simplifying assumptions 
can be used; constant thickness, constant angle (conical 
flows), or vertical equilibrium. Chakrabarti & Das (2001) 
proved that there is no essential difference among these 
three assumptions from the point of view of the existence 
of stationary solutions. The thickness of the flow is ap- 
proximated as a constant for the sake of simplicity, as in 
Nakayama (1992), Nobuta & Hanawa (1994), Blaes (1987) 
or MTK, although the approximate thickness H ^ Cs/Ok 
deduced from the balance of the vertical pressure force 
and the vertical gravity is not constant. 

The stationary flow is described by the conservation of 
mass and the Bernoulli equation: 



prvr 



const 



1 



2(r - 1) 



const 



(1) 



(2) 



where p is the density, Vr is the radial velocity, Cg is the 
sound speed, and I is the speciflc angular momentum. An 
example solution of unperturbed flow is shown in Fig. 1, 
for which I = 1.87 and Cg = 0.1. As a consequence of the 
inviscid hypothesis, stationary accretion flows onto a black 
hole are very sub-Keplerian in their outer parts. The ques- 
tion of the origin of this sub-Keplerian flow is still uncer- 
tain (Chakrabarti & Titarchuk 1995, Molteni et al. 2001). 
In Fig. 1, the unperturbed flow is supersonic between the 



outer sonic point (i?s2 = 44.5) and the shock (rgh = 6.8), 
becomes subsonic between the shock and the inner sonic 
point (i?g = 2.36) and goes into the central black hole 
supersonically . This example solution is precisely the 
one chosen by Nobuta and Hanawa (1994), in which they 
showed that the outer shock is stable to axisymmctric per- 
turbations. Our numerical results in Section 3, however, 
indicate that it is unstable to non-axisymmetric perturba- 
tions. 

The mass conservation equation and the Euler equa- 
tion are written as follows: 
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where w is the vorticity. In order to write the linearized 
equations in the simplest form, the two functions /, g are 
defined as follows: 

Sp 
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where / is the perturbation of the Bernoulli constant and 
g is the perturbation of the mass acc;rction rate. The fre- 
quency u>' of the perturbation measured in the rotating 
frame is defined as: 



w' = w — mfl , 



(7) 



where m is the azimuthal wave number and fl = l/r"^ is the 
angular velocity. With the standard method of linear sta- 
bility analysis, (perturbations proportional to 6""^*+*""^, 
see Appendix B), the following differential system is ob- 
tained: 
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where B = rVrWz, Wz is the vorticity along the rotation 
axis, M. = "Vr/cs is the radial Mach number, and the 
subscript 'sh' denotes the shock position. The boundary 
conditions corresponding to a perturbed shock velocity 
Awr are obtained in Appendix C: 
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Fig. 2. Eigenspectrum of the flow I — 1.87 and Cs — 0.1, show- 
ing the instability of the modes 1 < m < 5. The empty trian- 
gles correspond to the stable axisymmetric eigenmodes found 
by Nobuta and Hanawa (1994). 
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Fig. 3. The maximum growth rate of the non radial modes is 
reached for m 10 in the flow I = 1.87 and Cg = 0.1. The pat- 
tern frequency ujr/m of these unstable modes is concentrated 
in a limited band. 




where r] = ((ilogjV1/(ilogr)|sh (Eq. A.l), M and Vr are 

calculated at the post-shock side. 

The vorticity Swz produced by the shock is: 

5^^ = ^e^:.^'^''\ (13) 
rvr 

There are two differential equations for /, g and one un- 
known parameter lu in our system, thus three boundary 
conditions are needed to solve the equations. In addition 
to the two boundary conditions Eqs. (10-11) at the shock, 
a third equation is obtained from the critical condition at 
the sonic point, 

Wffson - w'/^on - ^^e-^'-sh =0. (14) 

These three boundary conditions are used to numerically 
solve the differential system Eqs. (8-9) and to determine 

the cigcnfrcqucncics lu. A single equation corresponding to 
this boundary value problem is formulated in Appendix D. 

3. Numerical results 

The standard Runge-Kutta method is used to integrate 
differential equations from the sonic point to the shock. 
Fig. 2 shows the eigenspectrum of the flow Z = 1.87 and 
Cg = 0.1 studied by Nobuta & Hanawa (1994), for pertur- 
bations < m < 5. The stability of axisymmetric pertur- 
bations confirms the results of Nobuta & Hanawa (1994). 
In contrast, non-axisymmetric perturbations are unstable. 
The highest growth rate is reached for m = 10 perturba- 
tions, as seen in Fig. 3. Determining numerically how the 
most unstable m depends on the value of {I, Cg) is beyond 




Fig. 4. The two thick solid lines are the threshold for shock- 
included solutions. The dotted lines measure the shock strength 

by the value of Aish indicated on the left. The thin solid lines 
correspond to the value of the minimum Mach rmmber A^min 
indicated on the right. The 19 flUed circles correspond to unsta/- 
ble outer shock solutions, and the 5 empty squares correspond 
to stable ones. 

the scope of this paper. In what follows, numerical calcu- 
lations are restricted to m = 1 perturbations for the sake 
of simplicity. 

Fig. 4 shows the domain (Z, Cg) of angular momentum and 
sound speed for which an isothermal inviscid flow, sub- 
sonic far from the accretor, may accrete onto a black- hole 
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Fig. 5. Growth rates compared to three timescales: advec- 
tion Tadv (squares), acoustic Tac (circles), and rotation f2(rson) 
(crosses). 



through a stationary shock. This domain is limited by two 
solid lines. The solid line on the left corresponds to the so- 
lutions in which the inner shock and the outer shock arc 
identical, whereas the solid line on the right corresponds 
to extremely weak shocks, A^sh ~ 1. A total of 24 shocked 
solutions were considered for numerical calculation, in or- 
der to get a broad view of their stability properties. Since 
the inner shock solution was already found unstable to 
axisymmetric perturbations, we concentrated on the sta- 
bility of the outer shock solution. Solutions at the bottom 
right corner in Fig. 4, for which the advection timescale is 
much longer than the acoustic one (TMmin < 0.01), were 
avoided for computing time reasons. Among this sample of 
24 shocked flows, 19 were found to be unstable and 5 to be 
stable. The results shown in Fig. 4 suggest that the shock 
is generally unstable to non-axisymmetric perturbations, 
except for a very narrow region close to the upper right 
border. In other words, the shock might be stable only for 
Msh ~ 1. 

4. Instability mechanism 

4.1. Comparison with the advective timescale 

In an isothermal flow, two types of cycles may exist be- 
tween the sonic point and the shock, i.e., the purely acous- 
tic cycle and the vortical-acoustic cycle. From the results 
of Foglizzo (2002), we would expect the vortical-acoustic 
cycle to be particularly unstable for very strong shocks, 
A^sh ^ 1, with a growth time comparable to the advec- 
tion timescale Tadv^ 
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Fig. 6. The position of the corotation radius is indicated for 
19 unstable flows (circles) and 5 stable ones (squares). 



The numerical results seem to exclude an explanation 
based on a vortical-acoustic cycle: 

(i) the flow is generally unstable to non-axisymmetric 
perturbations, even for mild shocks: the range of shock 
strengths among the 19 unstable flows considered is very 
wide (0.06 < M^h < 0.94). 

(ii) the growth time can be much shorter than the ad- 
vection time. The growth time of the instability is plotted 
in Fig. 5, in units of the rotation frequency at the sonic 
radius Cl{r son), the advection time Tadv and also in units of 
the acoustic time Tac approximated by the following lower 
bound: 



1-M 



2 

min 



(16) 



The dispersion of the points in Fig. 5 is smallest when 
the growth rate is measured in units of the acoustic time, 
suggesting a purely acoustic cycle. If this is the case, what 
is the mechanism ? 

4.2. The corotation region in the PPI mechanism 

The corotation radius rco of the perturbation is defined 
by a; = mflo, where Oq = ^(''co): 



Im 



(17) 



According to Fig. 6, the corotation radius of the most un- 
stable modes is always located between the sonic point 
and the shock. Fig. 5 and Fig. 6 are hints in favour of 
the PPI. The mechanism of the PPI was formulated most 
simply by Goldreich & Narayan 1985 (hereafter GN85) in 
thin discs. The more subtle effect of a corotation reson- 
nance (Narayan, Goldreich & Goodman 1987, Kato 1987) 
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Fig. 7. The value of the ratio Cs/^Kr at the sonic point (dotted 
lines) and at the shock (full lines) is indicated on the lines. 

can be neglected in the present study, since the epicyclic 
frequency is zero (r? = Z/r^). The PPI is based on the ex- 
change of energy and angular momentum between acous- 
tic waves propagating inside and those propagating out- 
side the corotation radius (Mark 1976). The thin disc hy- 
pothesis in those studies is technically important in order 
to treat high frequency perturbations in the WKB approx- 
imation uj = mOp ^ Cg/r. According to Fig. 7, the high 
frequency approximation is applicable only in the limit 
Cg — » 0, which coincides with strong shocks according to 
Fig. 4. Following GN85 and Kato (1987), we believe that 
the physical mechanism captured analytically in thin discs 
is also relevant in thicker discs {p-vsjca ~ 1), even if a 
quantitative estimate of the growth rate is precluded. 

The PPI occurs with either an inner or an outer re- 
flecting boundary, or even more efficiently with both. Let 
us recall this mechanism in the case of a single reflect- 
ing boundary, noting "T^o the fractional amplitude of the 
transmitted wave through the corotation zone, and the 
fractional amplitude of the acoustic wave reflected by the 
boundary. An incident wave with energy -1-1 transmits an 
energy — |7^oP on the other side of corotation, thus ampli- 
fying the energy of the reflected wave by a factor l-|-|7^oP- 
If T is the duration of the acoustic cycle, the growth rate 

of this instability is: 

1 



■log 



1^1(1 + l^coP)^ 



(18) 



The following estimate of T^o was obtained by GN85 and 
conflrmed by Kato (1987): 



T 



exp 



TT 
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m|(7o| 







(19) 



The instability is thus slow in thin discs if k > 0, and 
the most unstable modes correspond to m ~ Kr^ojcs- Let 



us now investigate the effects of radial advection on this 
mechanism. 



4.3. Effect of advection on the corotation region 

The effect of advection on the corotation region can be 
handled analytically in the high frequency limit. The sec- 
ond order differential equation satisfied by the homoge- 
neous solution /o is obtained in Appendix E: 



1 - 
M 



/oc 



-r 



, (20) 



In the high frequency limit, acoustic waves can propagate 
in the region where the function k'^{r) is positive, and are 
evanescent in the region where k'^{r) < 0. fc^ is approxi- 
mated at high frequency as follows: 



1 



(1-M2) 



(21) 



fc2 is positive near the sonic point {A4 ~ 1), negative 
near the corotation radius (w' ^ 0), and positive again 
at large radius. A calculation with k = and > in 
Appendix F shows that T^o is independent of the Mach 
number in the high frequency approximation, if the gra- 
dient of the Mach number is neglected: 



' exp 



2|r!o| 



(22) 



where CIq is the derivative of the rotation frequency at the 
corotation radius. This estimate at high frequency thus 
favours modes with a low azimuthal number m. Beyond 
the high frequency approximation, this calculation shows 
that the amplification mechanism at corotation can oper- 
ate even with radial advection. 



4.4. Crucial role of the boundaries 

In order to estimate the effect of radial advection at the 
inner boundary, Blaes (1987) considered an inviscid fiow 
with constant thickness in a pseudo-Newtonian potential, 
and a uniform angular momentum. It differed from the 
present study by an adiabatic hypothesis and, most im- 
portantly, a leaking outer boundary. As stressed by GN85, 
the PPI requires at least one reflecting boundary. The sta- 
bilizing effect of advection found by Blaes (1987) is thus 
directly related to the impossibility, for acoustic waves 
moving away from the corotation, to be reflected towards 
it. By contrast, a stationary shock standing at the outer 
boundary of the flow is a good reflector of acoustic waves, 
if the shock is not too weak. A continuity argument would 
lead to expect the reflection coefficient T^sh to decrease to 
zero in the limit A4sh ^ 1- This is confirmed by a calcula- 
tion of TZgh in the high frequency limit (see Appendix E) , 
which gives the same result as in isothermal flows without 
rotation (Fogfizzo 2002): 

1 - Msh 



(23) 
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Fig. 8. Growth rate measured in units of the total acoustic 
time Tac (circles) compared to the acoustic time Tout between 
the corotation and the shock (stars). 



According to Eq. (18) and (23), a strong isothermal shock 
(i.e. |7^sh| ~ 1) is a sufficient condition for the instabil- 
ity of the acoustic cycle between the shock and corota- 
tion in thin discs. Acoustic energy in shocked flows is thus 
trapped between the shock and the corotation radius, and 
may leak inside through the sonic radius without damp- 
ing the instability. On the basis of the results obtained by 
Blaes (1987), a strong decrease of the growth rate should 
be expected in the weak shock limit. This trend is clear in 
Fig. 5, when measuring the growth rate, as Blaes (1987), 
in units of the rotation frequency at a fixed radius. What 
is more, all the stable modes correspond to weak shocks 
{Msh > 0.8) according to Fig. 6. 

Fig. 8 shows that the dispersion of the points is further 
decreased by measuring the growth rate in units of the 
acoustic timescale Tout between the corotation and the 
shock, approximated by: 



Tout 



1-Mi 



(24) 



The combination of Eqs. (18), (22), (23) and (24) gives 
an analytical estimate of the growth rate which is justified 
only in the high frequency limit, i.e. for Cg and strong 
shocks. A direct application of this estimate outside its 
range of validity is unable to explain the strong disper- 
sion of the points in Fig. 8 for weak shocks. As noted in 
Sect. 4.2 with Fig. 7, solutions with a weak shock coincides 
with those with the largest parameter (cg/fixr), i.e. the 
least adapted to a high frequency approximation. It should 
be noted that part of the dispersion of the growth rate in 
Fig. 8 may come from the contribution of the vorticity 
perturbations, when the advection time is comparable to 
the acoustic time: they trigger additionnal acoustic waves 
which may add or substract to the acoustic power present 
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Fig. 9. Incidence of advection on the PPI for weak siiocks 
(A^sh = 0.9). The growth rate of the instability is measured in 
units of the acoustic time Tout between the corotation and the 
shock. 



in the acoustic cycle. Could this effect be essentially stabi- 
lizing ? According to Fig. 9, for a given shock strength (i.e. 
A4sh = 0.9), the fiow seems to be stabilized by a strong 
advection. 



5. Comparison with adiabatic simulations 

Although non-linear processes might redistribute the en- 
ergy of the unstable modes, it could be tempting to iden- 
tify the frequency w^™, observed in the simulations, with 
the real part of the most unstable cigcnmodc. Based on 
five regular oscillation examples, MTK found that the os- 
cillation period P is always in the same range, 60 < P < 
280 despite the large variation in the rotation rate at the 
shock distance, 100 < 27r/Osh < 2000. We find that aU 
the five regtdar oscillation examples listed in MTK satisfy 
the condition, Qgh < i^r™ = 27r/P < ^^son, which is a re- 
quirement of the PPI mechanism. 

A quantitative comparison with the results of MTK is 
hampered by our choice of an isothermal flow rather than 
an adiabatic one, and also by our restriction to m = 1 
perturbations which are not necessarily the most un- 
stable ones. A global understanding of which azimuthal 
wavenumber corresponds to the most unstable perturba- 
tion has not been reached yet. The only isothermal exam- 
ple considered in Fig. 3 with m > 1 perturbations favours 
8 < m < 12, which seems higher than the low azimuthal 
wavenumbers present in the simulations of MTK. 
Despite these uncertainties, the present isothermal calcu- 
lation is simple enough to illustrate the possibility of a 
PPI in a shocked accretion flow, occuring on a shorter 
time scale than any possible advective-acoustic cycle. This 
instability mechanism should thus also apply to adiabatic 
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flows, including those simulated by MTK, if the shock is 
able to reflect acoustic waves. The calculation ofTZgh in the 
adiabatic hypothesis (Foglizzo & Tagger 2000, Eq. (C13)) 
would be unchanged by rotation and non radial perturba- 
tions, in the high frequency limit, for a strong shock: 



n.. 



2i 



sh 



■7 



22 -|- 'y2 ('y — 1)! 



(25) 



Comparing Eq. (25) with Eq. (23), the reflection coef- 
ficient of strong adiabatic shocks is significantly smaller 
than strong isothermal ones in the high frequency limit. 
The adiabatic hypothesis in the work of MTK may 
thus marginally influence the strength of the instability 
(|7esh| ~ 0.36 for 7 = 4/3). 

A refined comparison of the growth rate and oscillation 
frequency with the results of MTK would require a linear 
stability analysis in an adiabatic flow involving also m > 2 
perturbations. This will be the subject of a forthcoming 
paper. 

6. Conclusion 

The main results can be summarized as follows: 

1. Despite post-shock deceleration, the outer shock 
solution is generally unstable with respect to non- 

axisymmctric perturbations. 

2. The growth time is comparable to the acoustic time if 
the shock is not too weak. The flow might be stable only 
when A^sh is close to 1. 

3. A WKB analysis of thin discs with advection proves 
that the acoustic cycle between the corotation radius and 
the shock is unstable if the shock is strong, thus extending 
the mechanism of the Papaloizou-Pringle instability to 
shocked discs with advection. 

4. Extrapolating the results established in thin discs, 
the instability of thicker discs could be interpreted as 
a form of the Papaloizou-Pringle instability modified 
by advection. This interpretation is supported by the 
existence of a corotation radius within the disc, a growth 
time comparable to the acoustic time, and the role of the 
shock strength. The influence of advection in the specific 
case of weak shocks, however, is still unexplained (Fig. 9). 



Appendix A: Description of the unperturbed flow 

We consider a simple one-dimensional steady fiow under 
the assumption of a constant thickness. The derivative of 
the Mach number is obtained by derivating the stationary 

flow equations (1-2): 



V 



dlogM 



(A.1) 



rflogr 

where Ik is the Keplerian angular momentum, = 
/2{r — 1)^. The continuity of the flow at the sonic point 
requires 



2 „2 , ,2 _ 'son 



. 



(A.2) 



The (/, Cg) parameter space for the shocked accretion flows 
is shown in Fig. 4. 

Appendix B: Linearized equations for 
perturbations 

The mass conservation and Euler equations are linearized 
in order to obtain the following three equations: 

dq im ■ i^P 

or r p 



.df . . I 



c„— iuSvr 6Wz = , 

or r 



imct 



f — itoSVy, + Vr6Wz = . 



(B.l) 
(B.2) 
(B.3) 



The evolution of vorticity is deduced from the curl of Euler 
equation: 

dw 

— +w{V-v) + {v-V)w-{wV)v^O. (B.4) 

The vorticity can be integrated when linearized, leading to 
Eq. (13). The differential system of Eqs.(8-9) is obtained 
from the above equations. The perturbation Sv,p is related 
to / and Ssh (Eqs. (B.3,13)), 



1 / r i^dr 

5v^ = - — imc^f + Bshe-^'-sh 



lojr 



(B.5) 



Eqs. (8-9), (B.5) enable us to get {6vr, Sv^, 5p) from (/, g). 



One particularity of this instability should be pointed 
out: although the instability mechanism relies on acoustic 
waves, vortical waves are continuously generated by the 
perturbed shock and advccted to the accretor (Eqs. 12, 
13). These vortical waves (or the entropy /vorticity waves 
in the case of non isothermal flows) may play the role of 
the 'blobs' invoked frequently to model the time variability 
of X-ray binaries (e.g. references in Belloni et al. 2002). 
The observed properties of these blobs could therefore be 
used to test the validity of a model of shocked accretion. 
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Appendix C: Boundary conditions at the shock 

The isothermal shock conditions are written as follows: 

P-Vr- = P+Vr+ , (C.l) 

p-{v^^_+cl)=p+{v^^+ + cl) , 



(C.2) 
(C.3) 



where the subscripts "-" and "+" denote pre-shock and 
post-shock quantities, respectively. Let the shock be per- 
turbed by in the radial direction: 



(C.4) 



Gu & Foglizzo: Non-axisymmetric instabilities in shocked accretion flows 



The perturbed velocity Au^ of the shock and its angle A6' 
with the azimuthal direction are 



dgo, 
dr ' 



I JO son son ) 

ulCs j w or 



^' dfo , 



AO 



dt 

1 (9 AC 

r dip 



A^ = AVr . 

r Lur 



(C.5) 



(C.6) 



' son (^son H" 1) 



4(r,on - 1)3 



(D.5) 
(D.6) 



Considering the velocity component perpendicular to the 
shock, in the frame of the shock, the isothermal shock 
conditions become the following: 

dVr- 



With the method in 'Appendix C of Foglizzo (2002), the 
following single equation for the boundary value problem 
is obtained: 



+ 



to 



-AC AVr = 

or uj ' 



f dfo Of 

J ' ~5 ' JO 

or or 



|sh 



P+ 



dr 



A^ + Sp 



sh 



/ /oAe 



dr 



where the parameter A is: 



X I 1;^+ + -^AC - —Avr + 5vr ) (P.7) 
or UJ 



X = 



+ ^AC 
or 



+ ^^AC Aw, 

or to 



r'^C^Vr 



iluj' 



— im — 



21 dlog(ri;r) 
dr 



(D.7) 



(D.8 



= [P+ 



dr 



■AC + 6p 



dVr 



dr 



-AC 



AVr + SVr I + C: 



(P.8) 



+ 



dv, 



and the value of df /dr at the shock is deduced from Eqs. 
(8) and (10-12). 

Appendix E: Reflection at the shock 

Perturbations fah,9sh at the shock are decomposed on 
acoustic waves fo,go and advected perturbations fB,gB- 



^A^ + Vr-Ae 
dr 



dv, 



dr 



AC + Vr+A0 + 5vJ^C.%) 



/sh = /o + /o + /b , 



(E.l) 
(E.2) 



rpi, , /inio\ A A AC The second order differential equation deduced from the 

ihe boundary conditions Eqs. (10-12) are deduced from /t-> t-> \ • 



the above three equations. 

Appendix D: Boundary value problem 

With three boundary conditions Eqs. (10-11), (14), the dif- 
ferential system could in principle be directly integrated 
using the Runge-Kutta method. The singularity at the 
sonic point requires to integrate from the sonic point to- 
wards the shock. The boundary value problem is conve- 
niently reduced to a single equation involving the solution 
fo of the homogeneous system. The homogeneous differ- 
ential equations are expressed as follows: 

Jo +50 



differential system (D.1-D.2) is 



Q^2 \ dr \ 



M 



2M 



c, 1-A^2 

fo 

2^2 



\dfo 
J dr 



duj' 



dr J (l-7W2)c; 



. 



(E.3) 



This equation is written in Eq. (20) in a canonical form 
as in Kato (1987), with 



fc2 = 



(1 



1 9^ l-M'^ 
2d^ °^ M 



dfo ^ 

dr Vr{l - A^2) 

dgo _ iLo' 

'd^ ~ Vril-M"^) 



■fo - M^go 



2 2 

im ct 



(jjr^Vr 



fo 



(D.l) 

(D.2) 



1 a , i-M' 



(E.4) 
(E.5) 



In the high frequency limit, the WKB approximations of 

At the sonic point, the following boundary conditions are ,i „ • , j-i- „„ i , j— „„„„„4-:„ „„i„i;^ „ „„„ 

T . J the ingoing jg and outgoing jg acoustic solutions are: 



obtained: 



poison — /o|, 

UJ 



(D-3) /o± oc 



LOM 



uj'p, J 



)■ 



exp 



r ■ ,j 



=F A* dr 



(E.6) 



The derivation at the sonic point can be obtained from 
L'Hospital method. 



The definition of the acoustic flux used in Foglizzo (2002) 
can thus be extended to rotating flows as follows; 



dfo. 
dr ' 



1 



2ilm 



"/o|i 



(D.4) 



F± 



Ma poj' 



(E.7) 
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The variations of /q are dominated by the phase varia- 
tions in the WKB hmit: 



fi 

dr ^ Cs 1 - A^2^° ' 



(E.8) 
(E.9) 



The dominant contribution of the advected vorticity to 

/ and g, at high frequency, is deduced from the differen- 
tial system (8-9) where the derivative were replaced by a 
multiplication by iw' /v^- 



/b = 



nii'z — Ijj' iB, 



sh 



^/2j,2 _|_ j7^2^2 g2 

imBsh 



(E.IO) 
(E.ll) 



Using Eqs. (E.IO, E.ll) and the boundary conditions (10- 
11-12), the linear system (E.1-E.2) is solved in order to 

obtain f^: 



AvrMil-M^)iMTt^y 



2Vr IJl. 



±1-IJlM 



i-qcs 1 - M"^ 
uj'r II 



(E.12) 



The reflection coefficient T^sh is defined using Eq. (E.7), 
so that |7^shP is the fraction of the reflected acoustic flux 
reflected inward: 



_ fl 
fo 



1 + i^M ( M — n 



1 - iJiM \M + n 



^ _|_ irjCs 1 — M 
io'r Ad — fj. 

w'r M+IJ. 



(E.13) 



(E.14) 



In the high frequency limit, required for the validity of 
Eqs. (E.6) to (E.ll), /x ~ 1 and the reflection coefficient 
is given by Eq. (23). 

Appendix F: Approximation of the amplification 
at corotation 

Defining a new variable X as in Foglizzo (2001), 

dX _ M 
"57 " 1 ' 



(F.l) 



the differential equation of acoustic perturbations is sim- 
ply 



5X2 



. 



(F.2) 



This equation can be approximated as a parabolic cylin- 
der differential equation in the vicinity of the corotation, 
by linearizing the rotation frequency and neglecting the 
variation of the Mach number in this region: 

Ma 



X 



1 



Ml 



(r - Tco) 



(F.3) 



where the subscript '0' denotes the corotation radius. The 
effect of advection is obtained formally from GN85 with 
K = by replacing the soimd speed Cg by Cg (1 — A^q) 2 , and 
the azimuthal wavenumber TObym/(l — A^q)^. The frac- 
tional amplitude of the transmitted wave, deduced from 
Eqs. (5-8) of GN85 is thus independent of the Mach num- 
ber in this approximation: 



%o = exp 



TT mCs 



2\no\ 



(F.5) 



An equivalent way to obtain the same result is the WKB 
approximation of the tunelling obtained by integrating the 
variation of the radial wavenumber, as in Kato (1987). 
Denoting by hl, tol the two zeros of k deduced from 
Eq. (21), 



Tco = exp(-$LL) 



$T.T, = 



1 



kdr 



roL 



m 



,'2 



Til 
TT mCs 



(1 - 7W2)c2 



dr 



(1-7^2)3 



2 \a 



(F.6) 
(F.7) 

(F.8) 
(F.9) 
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